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ANNALS OF MATHEMATICS. 

Vol. IV. February, 1888. No. i. 

A SOLUTION OF KEPLER'S PROBLEM FOR PLANETARY ORBITS OF 

HIGH ECCENTRICITY. 

By Dr. H. A. Howe, Denver, Col. 

The solution of Kepler's Problem involves the transcendental equation 

M=E—esmE, (i) 

in which M and E are respectively the mean and eccentric anomalies, while e is 
the eccentricity. 

Let F(E)=^\\^E—M— sin {E — iW)] 

= ^ sin» (^ - ^) + ^ sin' {E - M) 

+ ^sin'(i?-J/)+. ... (2) 

From (i) and (2) we easily get Grunert's equation, 

i+g zo&\M sm{E — \M)-ir F {E) 
i—e—sm\Mco^\E—\M) — F\E)- ^^J 

Equation (3) gives 

tan (£- i^) = i±i tan W' (T ^ ,) cos 1^^^ ^-^)- ^4) 

Assume that tan (£' — \M) = ^-i-^ tan ^M. (5) 

Substituting in (4) according to (5) and reducing, we obtain 

sME'-E) = ^^^^^.2FiE) (6) 

= k.2F{E). (7) 

Hence for moderate eccentricities, E' — E is approximately equal to 2F{E). 

Let 2->i = E' — M—sm(E' — M). (8) 

Then E' — E is likewise nearly equal to 2;y, and E' — M — 27^ is nearly equal 
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to E — M. From (2) we obtain by differentiation 

dF{E) = [\sm\E— M) + ^ sin* {E—M)+ . . . ]cos{E—M)d{E—M). (g) 

Equation (9) shows that a small change in the value o( E — J/ produces a much 
smaller change in F (E) ; therefore, we write with but slight error 

2 F{E) = £' — M— 27j — sin {E' — M— 2ij). (10) 

From (8), omitting terms of the seventh and higher orders, we have 

ry = j^-L. sin' {E' — M) + ^% sin' {E' — M). (i i) 

Suppose that E' is changed by an amount equal to 27] ; i. e. that dE ' = 25y. 
Differentiation of (11) and reduction of the resulting equation, neglecting powers 
of sin (E' — Af) higher than the fifth, gives 

dri = ^sm\E' — M) cos {E' ~ M); (12) 

. • . y] — drj = ^^ sin^ {E' — M) + -i-^ sin' {E> ~ M) 

— ^L sin' {E' — M) cos {E' — M). (13) 

The last two terms of (13) have nearly the same magnitude, and opposite signs. 

We therefore write fj — dq =^ ^^^ s\r? {E' — M). (14) 

Now since (12) is obtained on the assumption that dE' = 2)y, 

■q — dri = \\E' — M— 27] — sin {E' — M— 2r])]. 

Hence by (10) and (14) we may write 

F{E) = j\sin^E'~M). (15) 

Equations (15) and (7) give 

sm(E' ~E) = ^ksm^E' — M). (16) 

Let E' — M— ^k sin' {E' — M) be an approximate value of £ — M; then 
2F\_E' — M—^k sin' {E' — M)'] will be the corresponding value of 2F{E), and 
by (7), k . 2F\E' — M— \k sin' {E' — M)^^ will be the corresponding value of sin 
iJE' — E). Differentiation of (7) gives 

cos {E' — E)d{E' — E) 

= k\_\sxr?{E-M)^\sm*{E^M)-Y . . .^cos{E~ M)d{E- M) 

= k\_\s\n''{E—M)^\s\xx*{E — M)^. . .]d{E—M) 

= kvd{E — M). (17) 

\{ E — Mhft increased hy E' — E — \k sin' {E' — M), and the correspond- 
ing increment of E' — ^be denoted by J (E' — E), we have 

d{E' — E) = kv IE' — E — \k sin' {E' — M)\ sec (£' — E). (18) 

But E' -E^A{E' -E) = k. 2F[E' — M— ^k sin' {E' — M)] 

+ 2F{E'-E). (19) 
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Substitute from (i8) in (19), and add and substract ^k sin' (£' — M) in the 
second member of the resulting equation ; reduction gives 

(£' — £)[i-[~ Jtv sec {£' — ^)] = [ i + ki' sec {£' — £)] . \k sin» {£' — M) 

+ k . 2£\_£' — M— \k sm\E' — M)] 
— |/§ sin' {£' — iW) + 2F{£' — £). (20) 

Division of (20) by the coefficient of £' — £, and subtraction of 
2F{£' — £) from each side, gives 

sin {£' — £) = \k sin' {£' — M) 

k{2F\_E' — M—^ksirr{£'~M)']—:^sm^{E'~ M)} 
+ \ + k'c sec {£' — E) 

kv sec {£' — £) . 2F{E' — £) 

I -^kvsec{E' —£) ^^'^> 

But we have, with sufficient accuracy, 

2F {£' — £) = \ sin' {£' — £) = tAo'^ sin" {£' — M), 

kv sec {£' — £) . 2F {E' — E) _ k" sin"(^^' — M) 
^"•^ \-^kvs.&c{E' — £) "' ~ 2592ti +"fo) ■ ^^'^> 

We next determine the error of (21) when its last term is neglected and the 
quantity sec {£' — £) is dropped from the denominator of the second fraction. 
The eccentricities of the orbits of all but 19 of the first 268 asteroids are below sin 
15°. Of these 19, 15 lie between sin 15° and sin 20°, and 4 are greater than sin 
20°; the orbit of Aethra (132) has the greatest eccentricity, e being equal to sin 
22°. 5. On the assumption that e = sin 23°, we find that the maximum value of 
the last term of (21), as computed by (22), is o'.'oo3 -j-. If the factor sec {£' — E) 
in the denominator of the second fraction of (21) be dropped, an error of o'.'oo02 
may result. To compute 2' in the same denominator, we substitute for E — M 
(in the expression given for v in (17)) the value E' — M — ^k sin' {£' — M)_ 
Fromequations(i7)-(2i) wesee that, had E' — £ — \ksir?{E' — M), the incre- 
ment oi £ — M, been infinitesimal, E' — M — \k sin' {£ ' — M) might have been 
used for £ — Mm computing v, without error. But since E' — E — \k sin' 
(E' — M) is a small finite quantity, the proper value oi £ — M to employ is 

£'- M~\k sin' {£'- M)-\iE'-E~\k sin' {£' - M)]. 

The greatest error m E' — £ caused by using E' — M — \k sin' {£' — M) 
for £ — M, in finding v, is less than o'.'oo2. Greater accuracy is unnecessary, but 
could be attained, since i + k'v is always clo.se to unity, by using 

E< — M— \k sin' {E' — M) 

— \k {2F\E' — M— \k sin' {£' — MJ] — ^ sin' {£' — M)}, 
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or E' — M—^k sin' {E' — M) — kF \_E' — M— \k sin« {E' — M)\ 
Equation (21) may therefore be written 

sin {E' — E)=\k sin« (E' — M) 

k{2F{E'-M-\k sin» {E' -J /) ] - \ sin' {E' - M)} 

'^ " I + /fe • ^ ^^ 

When e = sin 23°, the maximum numerical value of the last fraction of (23) 
is 38'.'6. A simple approximate value of this fraction in terms of i and sin 
(E' — 3f) may be obtained as follows : — 

From (9) and (2), by putting E' — Mior E — M, ^k sin« {E' — M) for 
d{E — M), and I — i sin^ {E' — M) for cos {E' — M), we get 

2F[E' — M—^k sin' {E' — M)] = | sin' {E' — M) -\- (-^ — ^^k) sin'(£'' — M) 

+ (TfY-A^)sin'(£'-^). (24) 

Since kv is small, we assume that 

^-±^=,^kv=i-^ks\n^{E'-M). (25) 

Substitution of (24) and (25) in the last fraction of (23) gives for its value 
^{gk-io^)sin'{E' -M)^{-^{iSk-7k^)-^^{g^-io/^)-]s\r{'{E' -M) 

Since sin {E' — M) = e sin E' , neglecting powers higher than the fifth, the 
preceding expression becomes xir (9^ — \ot^^ f svc? E' . Because this involves 
the fifth power of e, its value is small for moderate eccentricities ; nevertheless, 
the expression is not sufficiently accurate to be employed in seven-place computa- 
tions with large eccentricities. As the second member of (23) is not easy to 
compute, it may be simplified, and its errors tabulated. It is not difficult to show 
that the maxinium value of k \%(\ — ^)~', which is reached when v ^o. The 
maximum of kv is near ^/, which is only 0.08 when e ^ svn 2-^° ; kv being, there- 
fore, always small, neglect of the denominator of the last term of (23) introduces 
into E' — .£ an error the value of which is not far from 

y^ kv {gk — lO/^^) ^ sin' E', or -^ {g^ — lO/^) / sin' E'. 

When e = sin 23° this error may arnount to 3"; when e = sin 15° it does not 
reach o'.'2. Rejecting i -\- kv and writing E' — E for sin {E' — E), (23) becomes 

E' — E =^ k . 2F\_E' — M — ^k sin' {E' — M)] (26) 

To expedite computation, two tables are needed, one giving the correction 
of the second member of (26), the other containing 2E{x) with the argument x. 
The computer having e and M given, would use (5) and (26), taking into account 
the tabulated correction to the value o{ E' — E given by (26). From the value 
of E thus found, the true anomaly could be obtained by one of the usual methods. 



